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We present an i?-parity conserving model of sneutrino dark matter within a Higgs-philic (7(1)' 
extension of the minimal supersymmetric standard model. In this theory, the ^ parameter and light 
Dirac neutrino masses are generated naturally upon the breaking of the 17(1)' gauge symmetry. One 
of the right-handed senutrinos is the LSP. The leptonic and hadronic decays of another sneutrino, 
taken to be the next-to-lightest superpartner, allow for a natural fit to the recent results reported 
by the PAMELA experiment. We perform a detailed calculation of the dark matter relic density in 
this scenario, and show that the model is consistent with the ATIC and FERMI-LAT experiments. 

PACS numbers: 95.30.Cq,12.60.Cn,12.60. Jv,13.85.Tp 



I. INTRODUCTION AND MOTIVATION 

The question of the nature of dark matter (DM) , which 
forms nearly one-fourth of the total mass in the universe, 
is a pivotal question for cosmology, astrophysics, and par- 
ticle physics. Other than its relic density little is 
known about the structure of DM. Recently, there has 
been great excitement that new information about DM 
may be revealed by the results from the PAMELA satel- 
lite claiming an increase in positron fraction in cosmic 
rays with energies above 10 GeV. Others, including the 
PPB-BETS [3 and HESS ^ experiments, also claim en- 
hancements in electron/positron flux for energies above 
100 GeV. More recently, the FERMI collaboration re- 
leased a measurement of the e+e^ flux in the 20 GeV 
to 1 TeV range All of these experiments are con- 
sistent with a new primary source contributing to local 
electron/positron fluxes in the 10 — 1000 GeV range 0. 
These anomalies in cosmic ray fluxes may well originate 
from astrophysical sources such as pulsars or nearby su- 
pernova remnants [7|. However, they might also result 
from the presence of DM particles, providing a golden 
opportunity to learn more about the properties of DM. 

This intriguing possibility has led to a number of pro- 
posals that interpret the data as arising from dark mat- 
ter annihilation Hjwhich necessitates an 0(100) boost 
factor) or decays [9| (which must proceed much slower 
than the present value of the Hubble parameter). More- 
over, PAMELA, while reporting an excess in the elec- 
tron/positron flux, claims no excess in the proton/anti- 
proton flux (lo| . Thus, one intriguing idea is that the 
DM carries lepton number, e.g., it is composed of right- 
handed scalar neutrinos (sneutrinos) in minimal |ll| - 
[isj and extended [l4| - [l6j supersymmetric models. This 
framework might explain the preference for leptonic an- 
nihflation/decay products [l7H23j . 

In the simplest extensions of the minimal supersym- 
metric model (MSSM) that include right-handed neutri- 



nos, the right-handed neutrinos are gauge singlets. Hence 
the strength of their Yukawa interactions with the stan- 
dard model (SM) flelds governs whether they thermalize 
together with the rest of the MSSM matter [H, [11 [11]. 
On the other hand, if the lightest sneutrino i/^ is the 
lightest supersymmetric particle (LSP) and therefore is 
a potentially viable DM candidate, it can be produced 
with the right relic density from the decays of heavier 
superpartners [ll|. If the heavier right-handed scalar 
neutrino vj^ is the next-to-lightest superpartner (NLSP), 
its lifetime turns out to be longer than the age of the 
universe [III [11] . Therefore, decays of vj^ should still be 
active, resulting in observable effects like the excess en- 
ergetic electrons/positrons in cosmic rays. In the MSSM 
with right-handed neutrinos, these features have been an- 
alyzed with respect to the observed cosmic ray anomalies 
[23| . However, such right-handed sneutrino dark matter 
scenarios generically suffer from two aesthetic problems: 
(i) the origin of the suppressed Dirac neutrino Yukawa 
couplings (though such small cou plings a re allowed on 
the basis of technical naturalness [llUl3l. [23j). and (ii) 
the origin of the supersymmetric /i parameter (23 |. In 
this paper, we present a ?7(l)'-extended MSSM model 
in which these aesthetic problems are remedied, and the 
desirable features of sneutrino DM, including the abil- 
ity to account for the PAMELA positron excess, can be 
correctly produced. 



II. THE HIGGSPHILIC (7(1)'-EXTENDED MSSM 

It is well known that the MSSM suffers from a natural- 
ness problem due to the presence of the superpotential 
bilinear operator • Hd [llj . Though the mass param- 
eter /i enters from the superpotential, it must be of the 
order of the mass terms associated with the supersym- 
metry breaking sector. This puzzle can be remedied by 
extending the matter and gauge structure of the MSSM, 
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e.g. within unified and/or string models [25| . To this end, 
theories with an extra U{1)' broken at the electroweak- 
to- TeV scale by SM singlets are known to be able to gen- 
erate an appropriately sized fi parameter (see e.g. [26|). 

The U{\)' symmetry can also play a crucial role in 
generating neutrino masses. The right-handed neutrino 
sector and the fx parameter can be correlated for both 
Majorana ^3] and Dirac masses 28 1. We assume here 
that lepton number is an accidental symmetry that is 
conserved at the perturbative level. Hence, the neutri- 
nos are Dirac fermions, requiring Yukawa couplings of 
O (IQ-^'^). These couplings are technically natural, but 
an explanation for such a strong suppression is clearly de- 
sirable. One way this can occur is if the 17(1)' invariance 
suppresses leading order contributions to Dirac neutrino 
masses and allows higher-dimensional operators [28 1. 



In this work, we assume that the U{1)' charges satisfy 
Q'h ^^Hd 7^ to forbid the bare /i term, and +Q'l+ 
Q'n 7^ to forbid a bare neutrino Yukawa coupling. After 
including an SM-singlet chiral superfield the relevant 
part of the superpotcntial takes the form [2^, [2^ 
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in which the C/(l)' invariance requires Q'j^ +Q'h ~^Q's = 
and Q'h^ + Q'^ + Q'n + Q's ^ 0. In the above, Mr is 
a large mass scale, and and are the Yukawa cou- 
plings responsible for generating the /x parameter and 
neutrino masses. Upon the breaking of the electroweak 
and U{iy gauge symmetries, the effective low energy pa- 
rameters include the /x parameter 



and Dirac neutrino masses 
1 



(2) 



Mr 



(5)(i/0)Y,^Y4(i/°)/sin/3). (3) 



The effective neutrino Yukawa coupling Y,^ leads to neu- 
trino masses in good agreement with experiment: 



3 X 10" 



2.8 X 10-3 eV' 



1/2 



for |Y,| - 1, tan/3 = {H^)/{H^) - 1, and 

(S) ~ 3 TeV , Mr ~ Mgut = 10^^ GeV . 



(4) 



(5) 



Hence, in this framework the n problem is resolved and 
appropriately suppressed Dirac neutrino masses are gen- 
erated upon J7(l)' breaking. Parameter values were cho- 
sen to obtain neutrino masses which are in the right 
range. In the next sections, we give more details of the 
spectrum constrained by the experimental data. 



III. COSMIC RAY ANOMALIES 

To account for the anomalous cosmic ray fluxes, there 
arise several additional constraints on the model. First, 



we assume the right-handed neutrino superfields are to- 
tal gauge singlets to avoid unacceptably fast decays that 
would otherwise occur via gaugino mediated processes. 
Second, the C/(l)' charges must be assigned such that 
the successful generation of the /i parameter and neu- 
trino masses are maintained. The option that we pursue 
in this work is to have a "Higgs-philic" U{iy, similar in 
spirit to [22j . 

Enforcing these constraints leads to the charge assign- 
ments displayed in TableU] Note that of the MSSM fields, 
only the right-handed up quarks U and the Higgs fields 
and S have nonvanishing C/(l)' charges. The charge 
assignment in Table |T] is anomalous as it stands. How- 
ever, all gauge and gravitational anomalies can be can- 
celled either by invoking family-dependent C/(l)' charges 
as in [29I (though in this case one needs to worry about 
constraints from fiavor violation), or by augmenting the 
matter content of the theory by sets of vector-like quark 
and lepton fields, as well as additional SM singlets that 
are charged under the U (1)' gauge symmetry. The details 
will be addressed somewhere else in our future studies. 

With this charge assignment, scattering processes that 
involve 77„, S and U arc influenced by the 11(1)' gauge 
boson (Z^) and the C/(l)' gaugino {Z') [s^. At hadron 
colliders, the right-handed up-type quarks (ur and cr) 
can undergo Drell-Yan annihilation through Z'^ to pro- 
duce Higgs fields. The decays of the up- type squarks also 
exhibit novel branchings due to Z' exchange [30| . 
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1/6 
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1/3 
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-1/2 
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1 


1 
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1 


1 


1 





Hu 


1 


2 


1/2 




Hd 


1 


2 


-1/2 





S 


1 


1 








TABLE L The quantum numbers of quark (Q,U,D), lepton 
(L,N,E), and Higgs (Hu, Hd, S) superSelds. The superpo- 
tential couphngs of quarks and charged leptons are kept as in 
the MSSM. 



We now explore the possibility that the excess positron 
flux observed in the cosmic ray data is due to the presence 
of sneutrino DM. To do so, we will assume that the soft 
supersymmetry breaking sector of the theory is such that 
the lightest right-handed sneutrino is the LSP and the 
next-to-lightest right-handed sneutrino Vr is the NLSP. 
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More explicitly, 



rriT-.i < m^2 < m- 



(6) 



where 



denotes the remaining superpartner masses, 



including the heaviest right-handed sneutrino, v^. The 
third heavy particle, the NNLSP, is chosen to be the light- 
est neutralino, Xi f-nd it requires special attention since 
it decays late enough to have some non-thermal contri- 
bution to the DM relic density. The hadronic Big Bang 
Nucleosynthesis (BBN) constraints will also be consid- 
ered for the late decaying particles in the model. We 
discuss these issues further in the next section. 



A. The DM relic density 

For the spectrum considered above, the DM has two 
components, the stable LSP and the practically stable 
NLSP which lives longer than the age of the universe. To- 
ether, their relic density should reproduce the WMAP 
ll| result. The right-handed sneutrino population forms 
during the cosmic evolution from the decays of the MSSM 
spectrum plus exotics, which, with a few exceptions, 
should already have reached thermal equilibrium. A 
given sparticle decays into right-handed sneutrinos with a 

rate F ^ |Yy| x (sparticle mass). Therefore, the conver- 
sion rate into right-handed sneutrinos is quadratic in the 
effective neutrino Yukawa coupling Y,y, and the growth 
of the LSP number density is fast enough to produce the 
observed DM relic density 

We expect the third sneutrino to also decay to the 
LSP at freeze-out and thus contribute to the relic density. 
However, in our chosen parameter space, this sneutrino is 
taken very heavy (quantified below), much heavier than 
other sparticles. It would likely mix more with the left- 
handed sneutrino and decay very quickly compared to 
the NLSP. Our detailed analysis shows that its life-time 
(~ 10~^ sec) is larger than its freeze-out time 10"^" 
sec) but its non-thermal contribution to the relic density 
is negligible since the branching ratio for decaying into 
LSP or NLSP is suppressed (less than 10~^^). 

The right-handed sneutrinos exhibit a non-thermal dis- 
tribution. To see this, recall that the right-handed sneu- 
trinos mix with the left-handed sneutrinos and there- 
fore interact with the gauge sector. For right-handed 
sneutrino annihilation, only a four-point interaction, a 
left-handed sneutrino, or a Higgsino exchange can con- 
tribute. The annihilation is out of equilibrium as long 
as the interaction rate is smaller than the expansion 
rate h ~ T^/Mpi. The rate of the four-point interac- 
tion F4 - |Y^|''r < h, since |Y^| - lO^^^^ xhe left- 
handed sneutrino and Higgsino exchanges also have rates 
Fp^, F^ < h, and hence the right-handed sneutrinos do 
not thermalize before the electroweak phase transition. 
This is maintained after the phase transition if the mixing 
with the active (left-handed) sneutrino is small jTll [sij. 

For a detailed calculation of the relic density, the model 
must be considered in detail and all the mixings, masses 



and branching ratios must be determined so that we 
would calculate the relic density by counting all possi- 
ble channels. Moreover, we need to evaluate the would- 
be relic density of the late decaying NNLSP which is 
computed with through conventional methods. We im- 
plement the model fully into the CalcHEP [s^l package 
program with the help of LanHEP [ssj program. Once the 
CalcHEP model files are provided, we use the MicrOMEGAs 
[13] software for calculating the NNLSP relic density. 

In the implementation, we consider a normal hierarchy 
for the neutrinos and chose the mixing parameters to be 
sin^ 26*12 = 0.87, sin^ 26*23 = 0.92, and sin^ 26*13 = 0.02. 
The parametrization of the mixing matrix is identical 
to the one in the quark sector and the additional CP- 
violating phases are taken zero. As the positron excess is 
proportional to the neutrino masses, we keep these non- 
zero, consistent with the following constraints: Amli — 



(7.59 ± 0.20: 

lo-^ev^ m 



X lO-^eV^ and 
. We take Mr, 



Am^2 = (2-43 ± 0.13) x 

/i, tan/3, and as free 



parameters and express (S) , Y^ and Y^ in terms of them. 

In this model the sneutrinos also mix, and the mixing 
matrix can in general be expressed as 
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rat- 


1'" ] 
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ml, / 

'^RR / 





(7) 



where i,j are the flavor indices and the matrix elements 
are given by 

+ {HI?Q'h^Q'l + {S)^Q'lQ's) 

+ {S?Q'^Q's) + {rn^jf+M% 
: (ml,, r 



mz 



In 



m; 




(8) 



Here Af^. and M^. are the soft mass terms and A^. are 
the trilinear couplings (assumed diagonal), and /i 
are given in ([3]) and ([2]), respectively. Note that due to 
the specific U{1)' charge assignments in Table U some of 
the terms in (jS]) are zero. The sneutrino mass eigenstates 
v\ 2 are given by 



1 - 



cos0~ sinO- ~ 

fL—UR. Vl—VR 

-sinO- ~ cos0~ ~ ] \vl 




with the left-right sneutrino mixing angles Oj 



1 / -I- TO~ii 

- arctan — ^^^^ ^ 

2 \ ra£i — rati 



(9) 



(10) 
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such that TO^i < m^i, i = 1,2,3. In the rest of the 
paper we will refer to the right-handed sneutrinos as 
^R' * ^ 1, 2, 3 for simplicity, but in the numerical analysis 
all the mixings are implemented. Note that even though 
such a mixing is small numerically, it must be retained. 
The alternative is to introduce mass insertions to obtain 
most of the decay channels discussed below and CalcHEP 
cannot handle such insertions. 

For the numerical analysis, we use the following input 
values: 

/i= 1560GeV, tan/3 = 0.15, M;^ = 10^^ GeV 
= 0.4, Ml, = 2000 GeV, A^, = 200 GeV 
Mat, =100GeV, Mat^ = 1.2TeV, Mn:, ^ 5 TeY 
Ml = 1400 GeV, M2 = 1600 GeV, M3 = 2TeV(ll) 

where Mi,M2 and are the soft gaugino mass terms 
for U{1), SU{2)l and SU{3), respectively. For the chosen 
point in the parameter space, the masses entering the 
calculations are given in Table [TTl Clearly, the LSP is ?g 
and the NLSP is vj^ , denoted as i/^ and t'^ in the rest of 
the text. The third right-handed sneutrino v'^ becomes 
ly^. The NNLSP is Xi ^.s previously chosen. The total 
widths of these particles are also given in Table IIIIl 

We first address the hadronic BBN constraints on sev- 
eral particles in the table, as it is known that late decay- 
ing particles may spoil the predictions of BBN and thus 
need special attention. From Table Hill there are only two 
particles, namely Xi and 1/^, with lifetimes ^1.2 sec and 
~ 10~^ sec, respectively, which survive and decay after 
the freeze-out of the NNLSP. In our scenario, the relevant 
quantity [s^ that needs to be evaluated and compared 
with the observational constraints on light elements like 
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m-o 
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m~o 
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m-o 
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1561 
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2537 


2538 


m~+ 
Xi 


m~+ 
X2 


m^fO 






niz' 


1513 


1648 


87 


2556 


5410 


2536 










m?^ 




1999 


1999 


1999 


1999 


1999 


1999 



TABLE II: The relevant masses in GeV for the chosen pa- 
rameter set. m^o (m.+ ) denotes the mass of the a^^ib*'') 

neutrahno (charguio) state. Hi 2,3 are the physical CP-even 
neutral Higgs bosons, m^i and m^i denote the masses of the 
LSP and NLSP refereed as m-i and m-2 ifl the text. 
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3.9 1.1 


X 10" 


i« 3.9 


4.6 X 10" 
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TABLE III: The relevant total decay widths for the particles 
in Table\l^in GeV for the chosen parameter set. F^i and T^i 

e fj. 

denote the decay widths of the LSP and NLSP. 



''He, D, ^Li etc is 

where e~i,2 are the energies carried by the hadrons and 
Y^o is the yield variable, defined as Y^o = n^o Is with n^o 

Xi ^ ' Xi Ai ' Xi 

the number density of Xi and s the total entropy density. 
This quantity is calculated with the help of MicrOMEGAs. 
We also generated all the relevant three-body decays 
eventually leading to hadronic final states. These are 
X? ^ and Xi ^ I'B^^t'^- Then H^, Z, and 

W decay hadronically. We found that B{H'l — > qq) = 
0.95, B{Z qq) = 0.72, and B{W qq') = 0.67. When 
all of the quantities are entered numerically, the variable 
■q is around 10~^^ GeV. In the range of the lifetimes of 
particles, the constraint coming from the overproduction 
of "^He is relevant. The observational constraint requires 
77(^He) < 10^^ '*^ [ll[, which is satisfied for the Xi case. 
We repeated the analysis for the heaviest right-handed 
sneutrino {v'^) and found r] even smaller than 10^^^ GeV. 

After insuring that the hadronic BBN constraints are 
satisfied, we proceed with the relic density calculation. 
As mentioned earlier, there are two particles, and 
as the components of the DM. Including the contribu- 
tions after freeze-out, the total relic density of the right- 
handed sneutrinos is given 

= nil + + ^'-°n (12) 

where "ce" ("fo") refers the contributions from decays 
at chemical equilibrium (after freeze-out). The contribu- 
tion to the relic density from the heaviest right-handed 
sneutrino is small and neglected here. Thus, £7^° has 
contributions only from the NNLSP Xa- It is defined as 

= ^S(x° ^ 1)], P.) + ^S(x l^l ^j) 

A 1 A 1 A 1 
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0.000112 


0.0444 


0.000929 


0.0665 


= 0.1119 



TABLE IV: The individual contributions to the relic density 
for the parameter set jllfl . 



Having shown that the parameter set of pT|) can re- 
produce the required rehc density, we next investigate 
the positron flux and discuss the PAMELA, ATIC and 
FERMLLAT data for the same parameter set. 



B. Understanding the PAMELA, ATIC and 
FERMI-LAT data 



where 1,2,3 (we also include the conjugated states). 
Of course, the sum of the branching ratios has to be unity, 
EL=i (^(x? + ^ ^,)) = 1. Nu- 

merically, we found Y^^i j=i^{Xi ~^ ^R^i) = 0-86 and 
Yllj=i^{x\ -> "i^RVi) = 0.14 in the parameter set ([TT|). 
fi^o is the would-be relic density of NNLSP for the case 
where it is stable. We have calculated Vl^o by using the 
MicrOMEGAs 34] through its CalcHEP interface. Numeri- 
cally, we found fij^o/i^ = 0.09 with the usual dimension- 
less parameter xp = Tp/m^o = 1/29.8. Here h and Tp 
are the Hubble parameter and the freeze-out tempera- 
ture. For example, the two most significant contribu- 
tions to from the (co-)annihilation channels are 

XiXi ^ tt (42%) and xtxi ^ tb (23%). 

For the calculation of ftl^ and , the following decay 

•"r 

channels are included (conjugated states are not listed): 



it- 



~2 



Viv'^^'^ , a=l,...6, 1 = 1,2,3 
' Iv]^' , 6=1,2, £^e,fi,T 

1.2 



W^^R 
^I'^R : 



= Bl, fJ,L,Tl,T2 
' ^I'^R ' 



(13) 



We performed a complete parameter scan, and found 
these are the only decay modes which are numerically 
significant. The individual contributions to the relic den- 
sity are summarized in Table IIVI In the table, we split 
the freeze-out contribution n~° into two parts, and 

r2p2 , which show the contributions of the NNLSP decay- 
ing to both the LSP and the NLSP , Xi i^i^w ^^^^ 
Xi — >■ f'ii^^) respectively. As can be seen from Table IIV[ 
the dominant contributions are coming from il-s (39.7%) 

'^R 

and il~°2 (59.4%) and only ~ 1% comes from decays to 

'^R 

the LSP, as expected. From ([TU]), the mixing angle is in- 
versely proportional to m~i — m% . Thus, as the mass 

'^LL _ ^RR 

difference between left and right handed sneutrinos be- 
comes smaller, their mixing becomes larger. From the 
fact that m^i = lOOGeV and m^2 = 1200GeV as well 

'^R '^R 

as mjji^ = 2000GeV, the NLSP V\ mixes largely with 
left-handed sneutrino fields and becomes most likely to 
be produced through the decays in ([T3| . The total relic 
density is 0.1036, which lies in the 2cr of the WMAP value 
(with those from the Sloan Digital Sky Survey) [37j : 



DM' 



0.111 



hO.Oll 
-0.015 



(14) 



Unlike the decays of heavy superpartners, the decay 
of the NLSP into the LSP proceeds much more slowly 
as seen from Table IIIII (see also [23] ) . The reason is as 
follows. As gauge singlets, the sneutrinos do not couple 
to gauginos but only to Higgsinos, via 



1 



5 UlK - eLH+) + HI^lS - H+blS 



'^vVR, 



(15) 



in which the Yukawa interactions of the Higgsinos are 



h^S [HuHd - H^Hj 



(16) 



These interaction terms show that the exchanges of 
and S induce 

Vr^v]iV,Vj, (17) 
and the iJ+ exchange gives rise to the decay 



^R ^ ^R^tq- 

The rate of this dileptonic decay is given by 



(18) 



(2 



i^(Y!:Y.)jY.Yl) 



m;72 / TO^2 



Ca 

22 \^sin/3 
2 \ 



32 \ TO~+ 



(19) 



after summing over all three lepton generations. Inte- 
grating over the Dalitz density gives 



Qi{x,y) = x{y-l)[x{b-Qx) + y{bx-2)] 
~ 2y^\ogy- \2{y-x){x-l) 

'x-y 



X iy ^ X + xy — 3x)\og 



x-l 



(20) 



which appears in (1191) with the indicated arguments. The 
two chargino states xt, ('^^ — 1j 2) with masses also 

appear in the decay rate via iJ+ = CiPrXi + C2PRxt ^ 
in which 



C1C2 _ ri4'2 (sin /3 -I- r^2 cos /3) 



CI -CI l-r2 +r2^2Cos2/3 



(21) 
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where C2 = ^/l - Cf, rw2 = V2Mw/M2, r^2 = 
hf^{S)/M2, and M2 is the SU{2)l gaugino mass. 

Various analyses [9;] of PAMELA and other satellite 
data suggest that, rather model independently, the DM 
candidate must have a lifetime of ^ 10^^ sec. This is 



consistent with the dilcptonic decay rate (|T9|) in that for 
the parameter choice in one finds 



r,±,T = 6.8 X 10"^^ GeV, 



(22) 



leading to a ^ 10^^ sec lifetime. This arises from the 
small (S) /Mqut ratio that sets the neutrino mass scale. 
One notices that the value of the rate involves the sum of 
all three generations of charged leptons (i, j = 1, 2, 3) as 
well as all the conjugated modes. Note that there is no 
contribution to the positron production from the third 
right-handed sneutrino since it only lives 10~^ sec and 
thus decays much faster. 

In contrast to the case of the MSSM with right-handed 
neutrinos [ill [23l | , the /i term and the neutrino Yukawa 
couplings are induced dynamically in our model. The two 
models also differ in terms of the predicted correlation be- 
tween the neutrino masses ([3|) and the dileptonic decay 
rate ([T9| . In the MSSM, these quantities are directly 
correlated via the Yukawa couplings ^ 0{1Q^^^). 
In our model, the rate involves both {S)/Mcut and 
mp^/mg , and, in the heavy Higgsino regime where 
m~+ ~ ~ /i ^ M2,Mw, it scales as (mj;|/mj;i^)'*, 

not as {{S) / Mgut)"^ ■ Hence, in the large /i limit there is 
only a mild dependence on , and the "see-saw scale" 
that suppresses the neutrino masses is then independent 
of the scale that governs the NLSP decays. However, 
in this regime the DM lifetime is longer than 10^^ sec, 
and hence the heavy Higgsino regime is not particularly 
preferred given present data. 

We now crystallize the qualitative approach presented 
above with a more detailed calculation of the positron 
flux originating from the decay (|T8l) with the decay width 
(ITOl) . We follow the procedure in to which we refer 
the reader for further details. The source term of the 
diffusion equation is 



Q^± {Ef.± , r) 



TNLSP dEf± 



(23) 



where tt-nlsp is the number density of the NLSP, tnlsp 
is the lifetime (which will be set as a free parameter to be 
fit), and dNg± /dEp± is the energy distribution of from 

the decay V'lf — vl^ifij. Before proceeding to evaluate 
dNg± /dEg± , we first determine tnlsp • 

Since both the LSP and the NLSP are components of 
the DM, other particles decay to either of them. Then, 
the ratio of the number densities of the LSP and NLSP 
should be equal to the ratios of the branching rates into 
either state. Explicitly, 



B'ul) 



(24) 



where A, B, A' , and B' represent all possible particles of 
the model and the right hand side is known once we fix 
the free parameters. For example, numerical evaluation 
yields, Y.^gB[A Bv]^) = 0.0718 and Y.a,b^{A ^ 
B vj^) = 0.428 for the particular parameter set considered 
here. There are two unknown in the left hand side of (|24p . 
The second relation needed arises from the mass density 
of the DM, pdm(''), as 



m^l "LSP + m^l '^NLSP = PDM{r) 



(25) 



where pDM(r) is further parametrized by adopting the 
Navarro-Frank- White density profile [s^. For the diffu- 
sion parameters, we use the three sets from models Ml, 
M2 and MED [|. 

We return to the evaluation of the dN^± /dE^± term. 
The definition and normalization of this term depends 
on whether the positron/electron pair is produced di- 
rectly from the decay — >■ v\^e~e'^ (that is, i = j = 1) 
or through some cascade decays. Our numerical anal- 
ysis shows that in addition to the direct e~e^ produc- 
tion case, the e^//^ has also significant decay width. 
In this case the muon /i^ will decay further as /^^ — 
i^/^e^i^e where the appropriate combinations of and 
Ue should be understood based on the charge of or 
e. In principle, the r lepton can also be produced di- 
rectly which then requires two cascade decays to get the 
positron/electron pair. So, we consider direct produc- 
tions of positron/electron, or indirect production only 
through the muon cascade. The details of how to treat 
the cascade productions of positron/electron from the 
muon cascade can be found for example in [39]. If 
the electron/positron pair is produced directly, then the 
dN^± /dE^± is given 



dNt 



1 



dE, 



dE, 



d^T 



(26) 



where Fg+g- is the decay width for the direct elec- 
tron/positron production case. If positron/electron is 
produced from the decay l/j^ — > i/^e^/i^ followed by the 
decay of the muon fi^ — >■ f^e^i^e, we write the dN/dE 
term for the positron only for notational simplicity (the 
electron case is very similar): 



dN 



dE,. 



dE„+ 



dE., 



dE^+ dE^j^+ 



27) 



where the second term is represented generically after 
carrying out the integral over the electron energy E^- . 
In this term, the muon decays further to a positron with 
energy _Ee+ . The combined contribution is included in 
the second term. The first term corresponds to the case 
where the leptons in the final state are e''"/i~ so that there 
is no further decay considered for the positron case. The 
normalization factor Ne±^^ is given by 



^2pe /.+ 



dE^+ dE^+ — — — — 



(28) 
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The solution of the diffusion equation for the flux $g± 
is given in terms of the Green's function in [9] as 



An rrieff tnlsp 



J dEG{E,E^±) 



dE 



(29) 



where either (gS]) or (H?]) should be used for the dN^± /dE 
term depending on the relevant decay mode. Here nries = 
a m^i + is an effective mass term originating from 
the solution of ((24)) and ([25t . We further approximate the 
Green's function as in [40] to maximize computational 
efficiency. 

The flux $g± obtained as a solution of the diffusion 
equation is then used to define the ratio for the positron 
flux R„± : 



<j,tot ■ 



(30) 



where <J>*°* is the sum of $e± and the background flux in 
our galaxy and we define <i>*°' 



(31) 



where represents the uncertanities in the e^(e^) back- 
ground. (k~,k+) = (0.7,0.9) is used in the numerical 
study. The following electron and positron backgrounds 
((j)P"'°^ Q^j. galaxy are approximated as in 

units of (GeVcm^ secsr)"^ 



^sec 

e+ 



0.16E- 



-1.1 



^ll£;o.9_^3 2i;2.i5 

o.7£;"-^ 

110£;i-5 + 600^2.9 
4.5^0-7 



580S4-2 



1 + 650£;: 



2.3 



1500£;4.2 



(32) 



The quantity i?e± in (1301) can be used to compare with 
the PAMELA data. Similarly the total flux <I>*°* + $*f 
scaled by is the relevant observable for the ATIC 
and the FERMI-LAT experiments. 

Our results for the positron flux are shown in Figs. ([1}, 
([2]), and ([3]) as the ratio i?e+ as a function of the positron 
energy Ef,+ for each propagation model Ml, M2, and 
MED, respectively. The direct and indirect production 
of the positron as well as the expected background are 
shown separately. The figures show that the models Ml 
and MED fit much better than the M2 model. The 
fitted curves significantly deviate from the PAMELA 
data (shown in the figures) for energies less than about 
5 GeV but are consistent with the experimental data 
for > 5 GeV. The fitted values of the lifetimes are 
0(10^^ sec), which are in the expected range. The life- 
times for the Ml and MED models are close to each other, 
and the value is only slightly smaller for the M2 model. 
Note also that the constraint on the lifetime of the decay- 
ing DM from the cosmic microwave background (CMB) 



Pi 



y 0.1 



0.01 



THLSP = 2.6 X lO^'sec 
Ml Model 



Backgrouncl 
PAMELA 




100 



1000 



Ee+ (GeV) 



FIG. 1: The positron flux as a function of positron energy 
in the Ml propagation model. The direct positron channel 
(e~e"'") and the indirect one through muon decay (fi^e^) are 
shown separately. The background and the PAMELA data 
are also shown. The fitted lifetime of NLSP is also indicated. 




0.01 



1000 



Ee+ (GeV) 



FIG. 2: The positron flux R^+ as a function of positron energy 
in the M2 propagation model. The direct positron channel 
{e~e~^) and the indirect one through muon decay (^^e^) are 
shown separately. The background and the PAMELA data 
are also shown. The fltted lifetime of NLSP is also indicated. 



is analysed in [43, and our results obey the lower bounds. 
The positron emitted through the /z^ channel deviates 
only from the direct channel part of the high energy tail. 

Simflar graphs are shown in Figs. dJ]), ([5]), and (|6]) 
for the total flux scaled by the positron energy E^± for 
the propagation models Ml, M2, and MED, respectively. 
Again the direct and indirect production of the positron 
as well as the expected background are shown separately. 
The ATIC and the FERMI-LAT data are also included. 
As seen from the flgures, our signal is a reasonably good 
fit to the ATIC and FERMI data. Our signal explains 
the FERMI-LAT data very well up to - 400 GeV for the 
e^e^ case and gets even better for the /i^e^ case, where 
the signal is consistent up to ~ 500 GeV. In the ATIC 
case, the signal explains the low energy as well as the 
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TNLSP = 2.4 X 10^'' sec 
MED Model 



^ Bacli ground 

e~e+ 

PAMELA 



1000 



0.1 - 




0.01 



100 



1000 



Ee+ (GeV) 



FIG. 3: The positron flux as a function of positron energy 
in the MED propagation model. The direct positron channel 
(e~e+) and the indirect one through muon decay (/x^e^) are 
shown separately. The background and the PAMELA data 
are also shown. The fitted lifetime of NLSP is also indicated. 




100 1000 

Ee+ (GeV) 



10000 



FIG. 4: The total flux scaled by E^± as a function of the 
positron energy in the Ml propagation model. Both the direct 
positron channel (e~e^) and the indirect one through muon 
decay (/i^e^) are shown separately. The background and 
the ATIC and FERMI-LAT data are also shown. The fitted 
lifetimes of NLSP for each data are also indicated. 



peak regions better than the intermediate enery interval. 
Once again, considering the fj,^e^ mode makes the sit- 
uation better in the end-region. In either case (ATIC 
or FERMI-LAT), the signal would fit better especially 
in the high energy region for another set of parameter 
values. For example, a bigger gap between the LSP and 
NLSP mass would allow more energetic which would 
shift the dying tail of the signal to the right. 

Unlike the PAMELA case, where Ml and MED sce- 
narios are favored, the M2 model works better for ATIC, 
while the FERMLLAT data slightly favors the Ml and 
MED models. The fitted hfetimes for the ATIC data are 
very close for the propagation models and are also quite 
consistent with the ones for the PAMELA case (espe- 
cially for the Ml and MED models). On the other hand, 
the lifetimes for the FERMI-LAT data are the same for 
the three models considered (the fitted curves are less 
sensitive to the types of the model used as long as the 
lifetimes are in the (5.4 — 6.0) x 10^^ sec range). For the 
Ml and MED models, the FERMI-LAT data requires 
lifetimes about twice bigger than than the ones for the 
ATIC and PAMELA cases. The ratio is bigger for the 
M2 propagation model. 

From the results, we can conclude that it is possible to 
satisfactorily explain the PAMELA-I-ATIC data simulta- 
neously (except the data in the 100-200 GeV range for 
the ATIC case) with a consistent lifetime for the DM. It 
is also possible to explain the FERMI-LAT data, which 
requires a slightly larger lifetime for the DM and works 
well for energies up to 500 GeV. 

Having discussed the electron/positron excess as well 
as the total flux, we now consider briefly the issue that 
no significant excess was observed in the anti-proton fiux 
[Toj . Unlike the case of the leptonic decays, the hadronic 
decays are generated only through the i^-term interac- 
tions in the superpotential. The relevant interaction 



terms in the Lagrangian are given by 



y!,y, 



21 \sin/3 



(^"')™^r^^K + H. c.,(33) 



1, 



in which {i 

bosons and TZui is the fraction of Hf in H^, [4 



, 3) are the neutral CP-even Higgs 

Higgs 



bosons produced in the vj^ decay [23| further de- 
cay into pairs of quarks (including top quark pairs as 
suggested by the mass spectrum in (fTT|) and in Table 
im gauge bosons, and Higgs bosons. The decay rate to 



1000 




100 1000 

Ee+ (GeV) 



10000 



FIG. 5: The total flux scaled by E''^± as a function of the 
positron energy in the M2 propagation model. Both the direct 
positron channel (e~e^) and the indirect one through muon 
decay (/i^e^) are shown separately. The background and 
the ATIC and FERMI-LAT data are also shown. The fitted 
lifetimes of NLSP for each data are also indicated. 
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10000 



Ee+ (GeV) 



FIG. 6: The total flux scaled by E^± as a function of the 
positron energy in the M2 propagation model. Both the direct 
positron channel (e~e^) and the indirect one through muon 
decay (/i*e^) are shown separately. The background and 
the ATIC and FERMI-LAT data are also shown. The fitted 
lifetimes of NLSP for each data are also indicated. 



quarks is 



gq 



(27r 



32 



Sq 



'M' 



H° 



,2 \ 



sin /? 



(34) 



The integrated Dalitz density has the functional form 



Gqix,y) = -4(1 -?/) + alogy + 



^(2/3^-7) 



X log 



7- (i-y)/3 
-y), = 



(35) 



xa + 7, and 7 



with a — —2a; + (1 
(l-y)2-a;(l + y). 

The light quarks are produced directly or indirectly 
through heavy quark, gauge boson, and Higgs boson de- 
cays, and hadronize to produce protons and anti-protons. 
An inspection of (|34|) and (fT9| shows that the anti- 
proton signal is suppressed compared to the positron sig- 
nal by a kinematic factor of order ^ 10'm'^/m~2 , which is 

O (lO"^*^) for u, d quarks and O (lO~^ ) fo r b quarks, as 
needed for consistency with PAMELA |l0| . 

However, the Higgs bosons that are produced in NLSP 
to LSP decays fragment efficiently into gauge bosons and 
top quarks, which in turn give off light quarks and anti- 
quarks with no apparent Yukawa suppression. Moreover, 
the leptons produced with the rate given in ([T9| can 
give rise to hadronic final states at the loop level. These 
indirect contributions are expected to yield anti-protons 
with an efficiency around 10% [H, HI]. 



The needed suppression of the anti-proton flux with 
respect to the positron fiux can stem from various effects, 
such as the interaction strengths in (|33|) . One way to 



(36) 



suppress Vqq with respect to Ti+i- is to have [22 
|(YtY.),J^<10-^ (YtY.)^^(Y.Yt)^^ 



which imposes an overall suppression on the fiavor- 
changing entries of the neutrino Yukawa matrices with 
respect to the flavor-conserving entries. A justification 
or realization of such a structure would in principle re- 
quire a detailed knowledge of the flavor structure of ([3|) . 
However, it is worth noting that there is no direct corre- 
spondence in general between the right-handed neutrino 
and the active (left-handed) neutrino sectors, which can 
lend credence for the needed suppression. 



IV. CONCLUSIONS 

We have presented a simple extension of the MSSM 
with an additional U{iy gauge symmetry that couples 
predominantly to Higgs fields. While other sneutrino 
dark matter models have been proposed to explain the 
cosmic ray observations, the model presented here does 
so while also resolving the naturalness problems of the 
MSSM. In this model, an electroweak scale fj, term and 
appropriately suppressed Dirac neutrino masses are gen- 
erated upon [/(!)' breaking. The right-handed sneutri- 
nos v]^ and Vj^ are the lightest and next-to-lightest super- 
partners, allowing for a natural explanation of PAMELA 
and other experiments like ATIC and FERMI-LAT in the 
context of sneutrino dark matter. 

A complete study of the DM relic density in this model 
includes contributions from the LSP and the NLSP, 
which is essentially stable with respect to the lifetime 
of the universe. We also analysed the contributions com- 
ing from late decaying particles, i.e., particles that decay 
after freeze-out. For the parameter set in ([TT]) . we obtain 
0.1036 for the total relic density of right-handed sneu- 
trino, consistent with the current WMAP value. 

We then discussed the PAMELA, ATIC and FERMI- 



LAT data, considering v\ 



^plflj as the source decay 
for the observed positron excess. The possibility of pro- 
ducing a final positron through muon is also discussed 
separately. The suppression of the anti-proton fiux com- 
pared to the positron fiux suggests the model-building 
constraint that the flavor-changing entries of the neutrino 
Yukawa coupling matrix are suppressed with respect to 
the flavor-conserving entries. The fltted lifetimes for the 
DM, of the order of lO^^ sec, flt the PAMELA data well 
for positron energies greater than 5 GeV, as well as the 
ATIC and FERMLLAT cases. The fit is better for the 
Ml and MED models described in the text. It can be 
concluded from the values of the fitted lifetimes that 
it maybe possible to explain the PAMELA-I-ATIC and 
FERMI-LAT data simultaneously in the corresponding 
energy ranges, though the FERMI-LAT data requires a 
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slightly larger lifetime values for the DM. In principle, it 
is foreseeable that a better scan of the model parameters 
as well as implementing an improved fitting procedure 
could offer an even closer simultaneous explanation of 
these three sets of experimental data. 
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